Optimally Designed Digitally-Doped Mn:GaAs 
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We use the ab initio local-density approximation (LDA) to study of exchange interactions and T c 
of Gai-zMn^As grown in digitally doped structures. We analyze the crystallographic dependence 
of exchange interactions predicted by the LDA in terms of the Mn t^ and e levels, and explain 
the origin of the antiferromagnetic contribution to the total exchange interactions. We exploit this 
dependence and consider 5-doping in specific orientations where the antiferromagnetic interactions 
are minimized, to optimize T c of the system. By including hole doping with the addition of Be in the 
GaAs host digitally doped Gai-^Mn^As is predicted to be significantly above room temperature. 

PACS numbers: 75.50.Pp, 75.30.Et, 71.15.Mb 
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Spintronics-based materials have attracted a lot of in- 
terest recently*»2. By introducing spin as a degree of free- 
dom existing technologies can be improved, e.g. smaller 
devices 3 - 

that consume less electricity, and are more powerful 
for certain types of computations as compared to present 
charge-based systems. A key point is that the magnetism 
be carrier-mediated, i.e. the magnetic state can be ma- 
nipulated by electrical or optical means. Dilute magnetic 
semiconductors (DMS), i.e. semiconductors doped with 
low concentrations of magnetic impurities Mn, Cr, or Co, 
especially Mn doped GaAs, are generally thought to be 
good candidates to meet the requirements for spintron- 
ics applications. A central issue that has impeded the 
spintronics field is the lack of a (demonstrably) carrier- 
mediated ferromagnetic semiconductor with critical tem- 
perature T c above room temperature. 

In Gai-zMn-rAs, Mn ions (approximately) randomly 
occupy the Ga sites and act as both spin providers 
(5=5/2) and acceptors. Because ferromagnetic cou- 
pling between local ions is mediated by carriers (holes in 
Gai-^Mn^As), the doping concentration is believed to 
be a crucial factor which determines the ferromagnetic 
properties of DMS. 

Traditionally, people use molecular beam epitaxy 
(MBE) methods to grow DMS thin films with 
magnetic elements randomly doped into the host 
semiconductors!**^. However, this preparation process 
always suffers of the low solubility of the magnetic el- 
ements, and the tendency of Mn to form interstitials 
deleterious to the magnetism. Although annealing sig- 
nificantly improves the quality of the samples, optimized 
Curie temperatures (T c ) for these random alloys are be- 
low 200K^, well below requirements for room tempera- 
ture applications. Although T c is predicted to be en- 
hanced when the material is co-doped with acceptors, 
attempts to do this with holes has not been successful. 7 
An alternative to a random alloy, <5-doping of Mn 
in a GaAs host has also been attempted with some 
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above the GaAs valence band, and is approximately lo- 
calized within the Mn layer, the holes are thus confined, 
which according to simple models should increase T c --*-. 

In this paper, we use a first-principles linear re- 
sponse approach based on the local-density approxi- 
mation (LDA), to study ferromagnetism in i5-doped 
Gai-^Mn^As, and contrast it to the random alloy. We 
show the following, and explain each point qualitatively 
from the electronic structure of Mn d levels: (i) The 
magnetic exchange interactions have a strong crystallo- 
graphic and concentration dependence. As a result: (ii) 
There is an optimal Mn concentration within the <5-layer; 
(iii) There is a optimal i5-layer thickness; (iv) p-type dop- 
ing (Be in GaAs) in the substrate can increase T c of the 
system by a large amount, bringing it well above room 
temperature. 

We adopt a linear-response approach!*! valid in the 
rigid-spin and long- wave approximation 13 , in which the 
change in energy from noncollincar spin configurations is 
mapped onto a Heisenberg form H = — J2rr' JrR' ' 
e.R>- 

In multiple-scattering theory the Jrr> are computed 
from the scaled susceptibility-!*-: 



Jrr< = — Im 

7T 



LL' 



doped in a GaAs epitaxial layer, produce high Mn doping 
concentration in a few (typically 2 to 3) GaAs monolayers 
(ML). Because the Mn-derived impurity band sits slightly 
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The notation and computational approach is based on 
the Atomic Spheres approximation and is described in 
detail in Ref. LLJ. Trl^li is related to the Green's func- 
tion and 8Prl is related to the exchange splitting. L 
is a compound index referring to the I and m quantum 
numbers of the orbitals at a site. 

Because of the large Mn local moment, the approxi- 
mations are expected to be reasonable. A previous pa- 
per using the same approach presents results for random 
alloys!*-. The finite temperature magnetization M (T) 
is obtained from the Heisenberg form assuming classical 
spins. The cluster variation method (CVM)!*- adapted to 
the heisenberg hamiltonian was used to obtain T c ; it was 
shown in Ref [15J that the (CVM) does a rather good job 
of estimating T c in random DMS alloys, and is expected 
to be comparably accurate here. 
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All the calculations are performed using supercells with 
periodic boundary conditions. The supcrcell we use in- 
cludes 8ML in the [001] direction. We set the thickness of 
Mn layer to be 1 or 2 ML, which is comparable with the 
experiment setup. The inter-layer distance is chosen to 
be 6 ML which is an optimized value based on both speed 
and accuracy considerations. Checks show that further 
increase of the interlayer distance does not significantly 
affect the exchange interactions T c . The lateral scale of 
the layers are chosen based on Mn concentration simu- 
lated. Both ordered and Special Quasi-Random Struc- 
tures (SQS)ii are used to simulate Mn doping inside the 
8 layer. 
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FIG. 1: Dependence of T c on Mn concentration inside an 
superlattice consisting of a 1-2 ML thick 5-layer of Mn, sand- 
wiched by a 6 ML thick GaAs layer. Mn are taken to form 
an ordered structure within the 5-layer. The x axis refers to 
the fraction of Ga atoms that are replaced by Mn in the layer 
considered. In the "m/n" label along x axis, m represents the 
number of ML in the Mn doped layer, and n is the total num- 
ber of cation sites in the doped layer. Circles show mean-field 
results for comparison. The 2/2 system is predicted to be a 
frustrated spin glass; this is reflected in the figure by assign- 
ing T c =0. Open diamond (circle) : T c for an SQS structure 
spread over 1ML at 50% concentration calculated within the 
CVM (MFA). 

We first investigate a series of ordered (5-layers varying 
concentration within the layer for 1-2 ML superlattices. 
Fig [1] shows the calculated T c for the Mn concentration 
in the 8- layer ranging between 1/4 and lJ£ 

For bulk random alloys^, the disorder in the Mn site 
positions rather strongly affects exchange interactions 
and reduces T c . In the bulk Gai-^Mn^As alloy, the pre- 
dicted optimum T c was found to be ^300K at £=12%. 
Fig. [T] compares an ordered to a disordered configuration 
of Mn within the <5-layer in the 1/2 filling case (Mn posi- 
tion in the ML were simulated by an SQS structure). It 
is seen that the effect of disorder is to reduce T c slightly, 
just as in the bulk random alloy. (Note that mean-field 
theory incorrectly predicts T c increases with disorder, as 
was also found in the bulk alloy.) 

Notably, the optimum T c for the 5-layer case occurs 
near 50% filling in a single ML and the optimum T c is 
near that predicted for the bulk case (~300K at x—12%). 



Most of this effect can be traced to the strong crystal- 
lographic dependence of the magnetic exchange interac- 
tions. This can be seen by considering exchange interac- 
tions Jrri in Gai-^Mn^As compounds of high Mn con- 
centration (Fig. [3]). While there is some variation be- 
tween compounds, there is a tendency for the Jrri to be 
weak or even antiferromagnetic along the [001] axis, but 
to be more strongly ferromagnetic along the [110] axis. A 
similar result is found for the dilute random alloys^. The 
origin of the strong crystallographic dependence can, at 
least in part, be explained in terms of the orbital charac- 
ter of the Mn d orbitals of t% and e symmetry. Mn impu- 
rities in GaAs form a "dangling bond hybrid" levels of 
t 2 symmetry which sits ~ O.leV above the valence band 
maximum, and is an admixture of the Mn d levels of t% 
symmetry (xy, yz, and xz) and effective-mass like states 
derived from the top of the host valence band. Below 
the valence band top sit Mn d orbitals of e symmetry 
(3z 2 — r 2 and x 2 — y 2 ), which hybridize less strongly with 
the environment than the t 2 orbitals. As is well known, 
the e states are filled, while the three ti states are par- 
tially filled with 2 electrons. (Mn, with 5 d electrons and 
2 sp electrons, substitutes for Ga, with 3 sp electrons. 
Thus Mn contributes 5 d electrons and a hole.) 
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FIG. 2: Dependence of T c on the Fermi level shift and excess 
hole concentration in a bulk SQS Gai^MnsAsioo alloy. Ef = 
corresponds to the usual charge-neutral case. 

While there is some disagreement as to whether the 
exchange is best described by an RKKY-like picture 2 ^ or 
a double exchange picture 21 , it is generally agreed that 
FM in Gai-^Mn^As derives from the partial filling of 
the t 2 level. According to the Zener-Anderson-Hasegawa 
mode l 22 ' 23 of double exchange/superexchange, exchange 
interactions between sites depend on the filling. At half 
filling the interactions are ferromagnetic, while filled lev- 
els couple antiferromagnetically. To the extent that this 
model is a valid description of FM in the Gai-^Mn^As 
system, the exchange interactions should be maximally 
FM when the t 2 level is half full. This is apparently a rea- 
sonable description of the LDA exchange interactions, as 
Fig. [U shows. Exchange parameters Jrri were computed 
in a bulk SQS-100 structure at 8% Mn concentration, as 
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FIG. 3: Scaled exchange interactions Jaw x J? 2 as a function 
of separation R = |R R' | . R is in units of the lattice constant 
a. (top) MnGa3As4 in the LI2 structure; (middle) the SQS-4 
Mn 2 Ga 2 As 4 ; (bottom) MnAs in the ZB phase. In the SQS- 
4 structure a given R be found in different environments, a 
rather strong dispersion in Jrw for a particular R is evident. 
The Table shows resolution of T C MFA into contributions from 
the Mn £2, e and sp orbitals. 



a function of Fermi level (thus altering the filling of the ti 
level), i.e. additional holes are modeled by simple Fermi 
level shifts. The Jrri and T c are maximal when the 
level contains approximately 3/2 electrons. 

If this applies to both the ti and e levels, the ti level 
should contribute a FM exchange, while the e level should 
contribute a competing AFM exchange. To see to what 
extent this model describes the LDA exchange we resolve 
the mean-field estimate for T c into separate / and m con- 
tributions, that is 



T, 



MFA 



MFA 
L 



T, 



MFA 



RL'R'L' 



(2) 



RR'L' 



where Jrlr'L' is identical to Eq. ([T]) but with the sum 



over LL 1 suppressed. Three compounds with high Mn 
concentration were studied: MnGa3As4 in the LI2 struc- 
ture, a four-atom SQS alloy Mn2Ga2As4, and MnAs in 
the ZB phase. Data shown in Fig. [3] largely verifies this 
picture. The table in Fig.[3jshows the e levels consistently 
contribute antiferromagnetically. FM contributions from 
ti level decrease with increasing Mn concentration. In 
the high concentration limit (ZB MnAs) the picture of 
a t2 derived impurity band is rather far removed from 
reality. Fig. [3] also shows site-resolved Jrri. It shows 
that Jrri when R points along [110] tend to be FM, 
while those pointing along [001] tend to be AFM. This 
tendency is compatible with symmetry of the ti orbitals, 
which point in the [110] directions and are FM coupled, 
while the e orbitals point along the [001] are weak or even 
AFM coupled. Thus it seems that the LDA exchange in- 
teractions are rather strongly dependent on symmetry of 
the Mn orbitals rather than largely being a function of 
the Fermi surface, as RKKY-like models suppose 20 . 

The orientation dependence of J has been studied 
previously by a number of authors, mostly in the di- 
lute case where the Mn doping concentration is < 
10% 24 i 25 i 26 i 27 i 28 . One approach estimates J from var- 
ious FM-AFM configuration a 24 ' 25 i 26 ' 29 . However, such 
calculations are predicated on the assumption that large- 
angle rotations can map to a Heisenberg model (a per- 
turbation theory valid for small rotations). More impor- 
tantly, this is approach is problematic especially for dilute 
case, because many J's are required. The longer-ranged 
Js are important to describe T c , where percolation is 
important^. The more rigorously based linear-response 
technique 27 explicitly calculates energy changes for small 
angle rotations. Many Js are calculated at once; and the 
range extends far beyond the unit cell. In any case, none 
of the preceding calculations yield AFM coupling in the 
[001] direction, because the competing antiferromagnetic 
interactions are much stronger at high concentrations 29 , 
at least in the LDA. In the dilute case, we obtain both FM 
and AFM exchange coupling for [001] direction, which 
when averaged yields a small net coupling, similar to a 
CPA linear-response calculation 2 ^. However, for high Mn 
doping in the (5-doped case studied here, AFM becomes 
increasingly dominant along [001]. 

The crystallographic dependence of Jrri , and the re- 
duction of Jrri with Mn concentration, help to explain 
why T c is maximal for partial fillings of the 5 layer 
(Fig. [1} , and also suggest that T c will be optimal for a 
rather small S layer thickness. To test this, we compute 
T c varying the number 8 layers. We fix the total thick- 
ness of the supercell as 8ML and change the thickness 
of (5-layer up to 8ML (corresponding to the bulk case). 
Each layer contains 50% Mn and 50% Ga (the 1-layer 
case corresponds to the "1/2" point in Fig. [T|). Data 
is shown in Fig. 2] T c is optimal at 2ML then decreases 
above that. (For thicknesses > 2 ML only T C MFA is shown 
CVM predicts a spin glass). 

We now turn to the question: can T c be increased by 
co-doping to add holes? We can mimic addition of holes 
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FIG. 4: Dependence of T c on the 5-layer thickness for an 
average Mn concentration of 50% per layer. Circles show 
mean-field results for comparison. 
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FIG. 5: Dependence of T c on excess hole concentration for a 
2ML (5-doped MnGaAs structure, with the S layer containing 
50% Mn, as calculated by the CVM. Hole doping is modeled 
through artificial Fermi level shifts, or induced through the 
addition of Be in GaAs layers adjacent to the 5 layer. 



by Fermi level shifts (as before); alternatively we can 
dope the GaAs layer with an acceptor such as Be. In 
Ref. HO, (5-doped MnGaAs were grown in this way. Extra 
hole's were supplied to the 5-layer by doping Be in the 
GaAs near (but not in) the 5-layer, to prevent compen- 
sating donors such as Mn interstitials from forming in 
the 5-layer. 



In Fig. El the dependence of T c on the number of extra 
holes is presented. T c increases markedly by the addi- 
tion of ~ 1/2 hole/Mn atom. Higher hole concentra- 
tions do not further increase T c (as in the random alloy). 
The difference between Fermi level shifts and using ac- 
tual dopants can be explained in that the latter contains 
electrostatic shifts as a result of charge transfer from the 
Be dopants to the Mn layers, thus reducing hole con- 
finement in the S layer. This confirms the experimental 
results for S layer and also suggest a much higher T c for 
defect free samples compared with the best T c recorded 
to date. 



In conclusion, we demonstrate a strong crystallo- 
graphic dependence of magnetic exchange interactions 
in Gai-zMurAs, and show that it can be exploited in 
(5-dopcd Mn:GaAs structures to optimize T c . We deter- 
mine the size and concentration of (5-layers, and hole con- 
centration that are predicted to result in an optimal T c . 
Under optimal conditions (which entails superlattices to 
be grown sufficiently defect-free) T c is predicted to be 
well above room temperature. 
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